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Contrast between Lagrangian and Eulerian analytic regularity 
properties of Euler equations 

Peter Constantin, Igor Kukavica, and Vlad Vicol 

Abstract. We consider the incompressible Euler equations on R'^ or T'^, where d £ {2, 3}. We 
prove that: 

(a) In Lagrangian coordinates the equations are locally well-posed in spaces with^erf real-analylicity 
radius (more generally, a fixed Gevrey-class radius). 

(b) In Lagrangian coordinates the equations are locally well-posed in highly anisotropic spaces, 
e.g. Gevrey-class regularity in the label ai and Sobolev regularity in the labels 02 , ■ ■ ■, aa- 

(c) In Eulerian coordinates both results (a) and (b) above are false. July 13, 2015. 


1. Introduction 

The Euler equations for ideal incompressible fluids have two formulations, the Eulerian and the 
Eagrangian one (apparently both due to Euler IIEul57ll ). In the Eulerian formulation the unknown 
functions are velocity and pressure, recorded at fixed locations in space. Their time evolution is 
determined by equating the rates of change of momenta to the forces applied, which in this case are 
just internal isotropic forces maintaining the incompressible character of the fluid. In the Eagrangian 
formulation the main unknowns are the particle paths, the trajectories followed by ideal particles 
labeled by their initial positions. The Eulerian and Eagrangian formulations are equivalent in a 
smooth regime in which the velocity is in the Holder class , where s > 1. The particle paths are 
just the characteristics associated to the Eulerian velocity fields. 

In recent years it was proved IIChe92L iGaiiiM ISer95L ISuelTl IGST121 IShnl21 lNadl3l iFZUl 
[ZFMICVWI41 that the Eagrangian paths are time-analytic, even in the case in which the Euler¬ 
ian velocities are only C®, with s > 1. In contrast, if we view the Eulerian solution as a func¬ 
tion of time with values in then this function is everywhere discontinuous for generic initial 
data lIC^ iHAMlOl IMY121 [MYUl This points to a remarkable difference between the Ea¬ 
grangian and Eulerian behaviors, in the not-too-smooth regime. 

In this paper we describe a simple but astonishing difference of behaviors in the analytic regime: 
The radius of analyticity is locally in time conserved in the Eagrangian formulation (Theorem ll.il) . 
but may deteriorate instantaneously in the Eulerian one (Remark 1 1.21) . Moreover, the Eagrangian 
formulation allows solvability in anisotropic classes, e.g. functions which have analyticity in one 
variable, but are not analytic in the others (Theorem 1 1.51) . In contrast, the Eulerian formulation is 
ill-posed in such functions spaces (Theorem 1 1.61) . 
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1.1. Velocity in Lagrangian coordinates. We consider the Cauchy problem for the incom¬ 


pressible homogeneous Euler equations 

ut + u-Vu+ 'Vp = 0 (1.1) 

V-u = 0 (1.2) 

u(x, 0) = uo{x) (1.3) 

where {x, t) e x [0, oo), and d £ {2, 3}. In order to state our main results, we first rewrite the 
Euler equations in Eagrangian coordinates. Define the particle flow map X by 

dtX{a,t) = u{X{a,t),t) (1.4) 

X{a,0)=a (1.5) 


where t > 0, and o € R'^ is fhe Eagrangian label. The Eagrangian velocity v and the pressure q are 
obtained by composing with X, i.e., 

v{a,t) = u{X{a,t),t) 
q{a,t) = p{X{a,t),t). 

The Eagrangian formulation of the Euler equations (ll.ll) - (ll.3l) is given in components by 

dtV^ + Y^dkq = 0, i = (1.6) 

Y^dkV^ = 0 (1.7) 

where we have used the summation convention on repeated indices. The derivatives dk are with 
respect to the label direction and Y^ represents the {k, i) entry of the matrix inverse of the 
Jacobian of the particle map, i.e., 

Y{a,t) = {VaX{a,t))-\ 

We henceforth drop the index a on gradients, as it will be clear from the context when the gradients 
are taken with respect to Eagrangian variables a or with respect to the Eulerian variable x. Erom 
(11.21) it follows that det(W) = 1, and thus, differentiating dtX = v with respect to labels, and 
inverting the resulting matrix, we obtain 

Yt = -Y{Xv)Y. (1.8) 

The closed system for (u, q, Y) is supplemented with the initial conditions 

v(a,0) = vo(a) = uo(a) 

Y{a,0) = I 

where I is the identity matrix. In the smooth category, the Eagrangian equations (ll.6l) - (ll.8l) are 
equivalent to the Eulerian ones (ll.ll) - (ll.3l) . 

1.2. Vorticity in Lagrangian coordinates. Eor d = 2 the Eulerian scalar vorticity uj = V"*- • u 
is conserved along particle trajectories, that is, the Eagrangian vorticity 

C{a,t) = uj{X{a,t),t) 


obeys 


C(o,f) = a;o(o) 


(1.9) 
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fort > 0. The Lagrangian velocity v may then be computed from the Lagrangian vorticity C, using 
the elliptic curl-div system 

SijYi^dkV^ = Y^dkV^ - Y^dkV^ = C = wq (1.10) 

Y^dkV^ = Y^dkV^ + Y^dkv"^ = 0 ( 1 . 11 ) 

where Eij is the sign of the permutation (1, 2) (*, j). The equation (11.101) above represents the 

conservation of the Lagrangian vorticity, while (11.111) stands for the Lagrangian divergence-free 
condition. Note that the right sides of (II.IOMI.III) are time independent. 

For d = 3 the Eulerian vorticity vector a; = V x u is not conserved along particle trajectories, 
and the replacement of (11.91 ) is the vorticity transport formula 

C\a,t) = dkX\a,t)4{a). (1.12) 

Thus, in three dimensions, the elliptic curl-div system becomes 

SijkYjdiv’^ = C = ( 1 - 13 ) 

Y’^dkV^ = 0 (1.14) 

where etjk denotes the standard antisymmetric tensor. In order to make use of the identity (11.131) . we 
need to reformulate it so that the right side is time-independent, in analogy to the two-dimensional 
case. Multiplying (11.131) with and summing in i, we get 


SijkYrYjdiv’^ = 


CO, 


0 > 


m = 1,2, 3, 


(1.15) 


which is a form of the Cauchy identity containing only Y. Recall here the standard Cauchy invari¬ 
ants IICau27llZFra 

EijkdjV^dkX^ = Wq, i = 1,2,3, (1.16) 

which can be obtained by taking the Lagrangian curl of the Weber formula IIWeb681 IConOll . Thus, 
for d = 3 we solve (11.141) and (11.151) for Vn in terms of Y and cuq. Note that, as in the d = 2 case, 
this system has a right side which is time independent. 

1.3. Isotropic and anisotropic Lagrangian Gevrey spaces. First we recall the definition of 
the Gevrey spaces. Fix r > d/2, so that is an algebra (we may replace with 

VF'’’P(M'^) for r > d/p and p G (1, oo)). For a Gevrey-index s > 1 and Gevrey-radius d > 0, we 
denote the isotropic Gevrey norm by 


m 


/3>0 m>0 


E 

=m 


l|9^/ll 


Hr 


(1.17) 


where /3 G Nq is a multi-index. Also, let Gg^s be the set of functions for which the above norm is 
finite. When s = 1 this set consists of analytic functions extendable analytically to the strip of radius 
6, and which are bounded uniformly in this strip (the latter property is encoded in the summability 
property of the norm). 

Similarly, given a coordinate 3 G d}, we define the anisotropic s-Gevrey norm with 

radius S > Oby 


m>0 


WdTfW 


Hr 


that is, among all multi-indices /3 with |/3| = m, we only consider /3 = {/3k) with (3k = mdjk, where 
6jk is as usual the Kronecker symbol. 
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1.4. Main results. We have the following statement asserting persistence of the Gevrey radius 
for solutions of the Lagrangian Euler equation. 

Theorem 1.1 (Persistence of the Lagrangian Gevrey radius). Assume that vq G and 

Vuq G Gs^s 

for some Gevrey-index s > 1 and a Gevrey-radius (i > 0. Then there exists T > 0 and a unique 
solution {v,Y) G C{[0,T]] x C([0, T], i^'') of the Lagrangian Euler system (ll.6l) - (ll.8l) . 

which moreover satisfies 

Vv,Y eL^i[0,T],Gs,s). 

On the other hand, if the uniform analyticity radius of the solution u{x, t) of (ll.ll) - (ll.3l) is 
measured with respect to the Eulerian coordinate x, then this radius is in general not conserved in 
time, as may be seen in the following example. 


Remark 1.2 (Decay of the Eulerian analyticity radius). We recall from Remark 1.3 in HKVllbll 
that there exist solutions to (ll.ll) - (ll.3l) whose Eulerian real-analyticity radius decays in time. Con¬ 
sider the explicit shear flow example (cf. IIDM87llBTIO]l i given by 


u{x,t) = {f{x 2 ), 0 ,g{xi -tf{x 2 ))) 


(1.18) 


which satisfies (11.1 Ml.21) with vanishing pressure in d = 3, for smooth / and g. Eor s 
for simplicity consider the domain to be the periodic box [0,27r]^, and let 


f{y) = sin(y) 


and 


9{y) 


1 

sinh^(l) + sin^( 2 /) 


1 we may 


It is easily verified that the uniform in xi and X 2 real-analyticity radius of u(x, t) decays as 


1 

f + 1 


for all t > 0, and is thus not conserved. Note however that the above example does not provide the 
necessary counterexample to Theorem ll.il since g does not belong to the periodic version of Gip. 
Indeed, (— 1)"'5^^"'^(0) > (2n!)/4, and thus the series defining ||p|| Gi 1 is not summable. 


The next statement shows indeed that Theorem 1 1.1 1 does not hold in the Eulerian setting. 


Theorem 1.3. There exists a smooth periodic divergence-free function uq such that 



||uo||gi,i < 00 

(1.19) 

and such that 




lk(i)llGi,i = 00 

(1.20) 

for any f > 0. 




The example proving Theorem 1 1.3 1 is provided in Section [5] 

We are indebted to A. Shnirelman nshuisn for the following example, pointing out that the 
results of Theorem 11.11 are sharp in the sense that the time of persistence T of the Eagrangian 
analyticity radius may be strictly less than the maximal time of existence of a real-analytic solution. 

Remark 1.4 (Time of persistence of the Lagrangian analyticity radius BShnlSIl ). Consider 
the stationary solution 

u{xi,X 2 ) = (sinxi cos X 2 , — cos xi sin 3 : 2 ) 
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of the Euler equations in This is an entire function of {xi,X 2 ), and moreover, the xi-axis is 
invariant under the induced dynamics. Abusing notation we denote by 

Xi(oi,f) = Xi(ai,0,f) 

the image of the point (oi, 0) under the flow map at the moment t, and by 

122 ( 01 , f) = 122 ( 01 ,0, t) = (9ajXi)(ai, 0, t) 


its Lagrangian tangential derivative. These functions satisfy the ODE 
^Xi{ai,t) = sin(Ai(ai,f)), Xi(ai,0) = ai, 
^122(01,0 = cos(Xi(ai,f))l 22 (oi,t), 122(01,0) = 1 . 


The solution Xi is given by 


cosXi(ai,f) = 


(e^* + 1) cos(ai) — (e^* — 1) 
(e^* + 1) — (e^* — 1) cos(ai) 


and its tangential gradient obeys 

122(01, f) = 


2e^ 

(e^* + 1) — (e^* — 1) cos(ai) 


Thus, for any fixed f > 0, the function 122 ( 01 , f) has a singularity at the complex point ai = 
3?ai + Sroi (and its conjugate) satisfying 


cos(ai) = cos(3?ai + i^ai) = 

so that 

3?ai = 0, and l^roil = In 
Note however that this singularity obeys 


e^* + l 

g2t _ I 

e* + 1 
e* — 1 


|ai| —>■ 00 as f 0^ 


( 1 . 21 ) 


In summary, at any fixed f > 0 the function Y (a, t) is not anymore entire with respect to the label 
a. Given any 5 > 0, we have XaVo = V^no € Gi,S< and while Xxu{-,t) = XxUq £ Gi^s for all 
f > 0, there exists 

^ = U«) = ln(|^) >0 
such that Y (•, f), and thus also Xav{-,t), obey 

l|l"(-,f)llGi,i, ||Vn(-,f)||Gi .5 ^ 00 as t^T{6)~. 

Thus, the time of analyticity radius persistence T guaranteed by Theorem 11.11 cannot be taken as 
infinite. Yet, Theorem 1 1.1 1 is consistent with T(5) —)• 0 as <5 —)• 00 . 

The proof of Theorem 1 1.1 1 may be used to obtain the local existence and the persistence of the 
radius for anisotropic Gevrey spaces as well. 
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Theorem 1.5 (Solvability in Lagrangian anisotropic Gevrey spaces). For a fixed direction 
j € {1,..., d}, assume that vq € and that 

for some index s > 1 and radius d > 0. Then there exists T > 0 and a unique solution 
{v,Y) G (^([0, T],F[^~^^) X (^([O, T],F[^) of the Lagrangian Euler system (ll.6l) - (ll.8l) . which more¬ 
over satisfies 

Vi;,FGL°°([0,r],Gg). 

The above theorem does not hold in the Eulerian coordinates as shown by the next result. The 
fact that the Eulerian version of the theorem does not hold might not surprise, due to the isotropy and 
time-reversibility of the Euler equations. On the other hand, the fact that the Eagrangian formulation 
keeps the memory of initial anisotropy is puzzling. 

Theorem 1.6 (Ill-posedness in Eulerian anisotropic real-analytic spaces). There exists T > 
0 and an initial datum uq G C'°°(M^)/or which uq and loq are real-analytic in xi, uniformly with 
respect to X 2 , such that the unique C([0, T]; H'') solution uj{t) of the Cauchy problem for the Euler 
equations (fTB-lfOl) is not real-analytic in xi, for any t G (0, T], 


2. Ill-posedness in Eulerian anisotropic real-analytic spaces 


In this section we prove Theorem IE6I Here, all the derivatives are taken with respect to the 
Eulerian variables. The idea of the proof is as follows. We consider an initial vorticity that is 
supported in a horizontal strip around the xi axis and which is nonzero in a horizontal strip and 
is very highly concentrated near the origin. We can construct it such that it is real analytic in xi, 
but is obviously not real analytic in X 2 - Given that the vorticity is approximately a point vortex at 
the origin, the corresponding velocity is approximately a pure rotation. Then for short time, the 
Euler equations will evolve in such manner that the vorticity is supported in a slightly deformed but 
rotated strip. The rotation uncovers some of the points that were on the boundary of the original 
strip, making them points of vanishing vorticity, while covering others. Thus, on a horizontal line 
parallel with the xi axis, the vorticity instantly acquires an interval on which it must vanish, while 
it is not identically zero, and hence it cannot possibly continue to be real analytic with respect to xi. 

In the detailed proof we first construct a function 

u{xi,X 2 ) G C°°(M2) 


such that the following properties hold: 

(i) div u = 0 on M^, curl u = oj, 

(ii) supp w C {(yi,2/2) : -1 < 2/2 < 1} 
(hi) U 2 (l, 1) > 0 and U 2 (—1,1) < 0 
(iv) There exists e G (0,1/2) such that 


uj{xi,X2) / 0 , (xi,X2) G {(2/1,2/2) : I2/1 - 1| < e, 1 - e < 2/2 < 1 } 

(v) (tangential analyticity for u) There exists constants Mq, 5o > 0 such that 

lom , M ^ Moml 

1^1 u{xi,X 2 )\ < 


with 


\&fdiu[xuX2)\ < 


Mom! 




m 

0 


(2.1) 


( 2 . 2 ) 
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and 


\d'^d2u{xi,X2)\ < 


Mqw}. 

"o 


(2.3) 


(vi) 5“ti; converges to 0 exponentially fast and uniformly as xi ±oo, uniformly in X 2 - 
In order to simplify the presentation, we introduce the following notation: If a; is a function (or 
a measure) with a sufficient decay at infinity, denote 


u{u;) = / K{x - y)uj{y) dy 


where 




denotes the Biot-Savart kernel. Now, choose a test function 


lA G Co°°(M) 

with values in [0,1] such that ^ ijj = 1 with > 0 for x G (—1,1) and ijj 
Consider the sequence of vorticities 

= cqA:^ exp(-A;^(xi + X2))iA(x2) 
for /c = 1,2 ,..., where cq is a normalizing constant such that 

/ oj^^\x) dx ^ 1 ask ^ oo 


Denote by 

u^^\xi,X2) = u{uj^^\xi,X2)), k = 1,2,... 


1 on [-1/4,1/4]. 

(2.4) 


the corresponding velocities. Each individual member of this sequence of velocities satisfies fhe 
assumptions (i), (ii), (iv), (v), and (vi). (Nofe however fhaf fhe consfanfs in (12.1I) - (I2.3I) depend on 
k.) The consfrucfion of a desired vorficify is complefe once we show fhaf for k large enough, we 
have 

4")(i,i) >0 

and 

4")(-i,i) <0. 


These inequalifies for k sufficienfly large indeed follow immediately once we observe fhaf fhe se¬ 
quence (12.41) is an approximafion of identify, i.e., if converges fo fhe Dirac mass (5o, while fhe velocity 


u = u{6o) 


J_ (\ 

271 y |xP’ |xP/ 


corresponding fo (5o satisfies (iii). Thus fhe consfrucfion of a velocity safisfying fhe properfies (i)- 
(vi) is complete. Denofe fhis velocity by uq and fhe corresponding vorficify loq = curlttQ. Now, 
consider fhe Euler equafion 


ujt + u(lo) ■ Vlo = 0 


wifh 


w(0) = Wo 
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where, recall, u{uj) denotes the velocity computed from the vorticity io via the Biot-Savart law. By 
the well-known properties of the Euler equation, the solution is smooth for all f > 0. By (ii) and 
(hi) and using the Lagrangian variables to solve the Euler equation, there exists fo > 0 with the 
following property: Eor every t £ (0, to), there exists a constant ei(f) > 0 such that 

uj{xi,X 2 ) = 0, |(xi,X2) - (-1,1)1 < ei(f) (2.5) 

On the other hand, by (hi) and (iv), we obtain, by possibly reducing to, that for every t £ (0, to) 
there exists a constant 62 (f) > 0 such that 

a;(xi,X2)/0, |(xi,X2) - (1,1)1 < e2(f)- (2.6) 

The properties (12.51) and (12.61) contradict the tangential analyticity of uj{t) at 0:2 = 1 for all t G 
( 0 ,fo)- 


3. Local solvability in Lagrangian anisotropic Gevrey spaces 


In this section we prove Theorem 1 1.5 1 Eor simplicity of the presentation, we give here the proof 
for d = 2. The proof carries over mutatis mutandis to d = 3, where the only change arises from 
using (11.151) instead of (11.101) . These details may be seen in Section IH where the well-posedness 
(by which we mean the existence and uniqueness) in 3d isotropic Gevrey spaces is proven. 

Eix s > 1. Without loss of generality, the direction j G {1, 2} may be taken to be j = 1. Eix 
(5 > 0 so that Vvq G with the norm M, that is, the quantity 

= WOT^voWnr 


obeys 


m>0 


V 12m—17 < M 
mx 


(3.1) 


Recall that Tq = I. 

Eix T > 0, to be chosen further below sufficiently small in terms of M, s, and 6. Eor m > 0 
we define 


Vm = Vm{T)= sup \\dY^Vv{t)\\Hr, (3.2) 

iG[0,T] 

= Z^{T) = sup t-^/^\\d^{Y{t) - I)\\Hr. (3.3) 

te[o,r] 

Observe thaf in the norm (13.21 ) the velocity v does not appear without a gradient. Also, we note that 
the power —1/2 of t appearing in (13.31) is arbitrary, in the sense that the proof works with any power 
in (-1,0). 

Eirst we bound Vu from the approximate curl-div system (ll.lOMl.l II) . in terms of Y and wq. 
Since (9™ commutes with curl and div, we may use the Helmholtz decomposition to estimate 


||(9™Vu||h'- < GIIOJ" curl + C\\d'^ Aiy v\\h^ . 

Eurther, by appealing to (II.IOMI.III) . the Eeibniz rule, and the fact that is an algebra, we obtain 

W^vWnr < C\\dT{u:o + eij{6^k - Y^)dkV^)\\Hr + C\\dT{{Sik - Y^)dkV^)\\Hr 

< C\\dTu,o\\Hr + C||y - I\\Hr\\dTVv\\Hr + C\\dT{Y - I)\\Hr\\Vv\\Hr 



\\di{Y - I)\\H4dT-^^v\\Hr. 
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Taking a supremum over t G [0, T] and using the notation (I3.2I) - (I3.3I) . we obtain 

Vm < CQm + + CT^/^ZmVo + ^ j ZjV,n-j (3.4) 

j=i V.?/ 

for all m G N, while for m = 0 we have 

Fo < cno + CT^/^ZoVo. (3.5) 

Note that we have not used here the evolution equation (11.61) for v, and have instead appealed to the 
Lagrangian vorticity conservation (11.101) . 

In order to estimate Z^, we use the Lagrangian evolution (11.81) in integrated form, and obtain 


I- 


Y{t)= /V: 

40 


Vn: Y dr 


[ {Y -I): Vv: {Y - I) dr + [ {Y - I): Vv dr 
Jo Jo 


+ 


[ Yv: {Y -I) 
Jo 


(3.6) 


— I) dr + / Vn dr 

lo 

for all t G [0,T]. Dividing by and taking a supremum over t G [0,T] it immediately follows 
from (13.61) that 

Zo < + T^^^ZofVo. (3.7) 

Differentiating (13.61) m times with respect to the label oi, using the Leibniz rule, and the fact that 
is an algebra, we arrive at 

\\dT{Y{t)-I)\\Hr 

rt 


< 


m 
j k 




^ dr 


E I 

+ 2 E / ( 7 ) " I)\\H4dT~'^v\\Hr dT + WdTYvWnr dr 

for all m > 1. Further, dividing by taking a supremum over t G [0, T] and using the notation 

(I3.2I) - (I3.3I) . we obtain 


\{j,k)\<m 


m 

j k 


j=0 


ZjZkVm-j-k + ( J ) i 


< CT^/^TZlVm + TZmZoVo + T^^^ZoV^ + T^^^Z^Vq + K^) 


m 
j k 


m—1 


ZjZkVm-j-k + CTj2i^.) ZjVr 


+ 

0<\ij,k)\<m 

for some constant C > 0. 

From (13.51) and (13.71) we obtain that for any t G (0, T] we have 
Vo(f) < CoDo + Zo{t)Vo{t) 
Zo{t)<Co4/‘^ sup (vb(r)(l+ f^/^Zo(T))^') 

TG[ 0 ,i) ^ ^ 


m-j 


(3.8) 
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for some constant Cq > 0 , while the initial data obey 

^o(O) = \\Vvo\\m=no<M 

Zo(0) = 0. 


Here we used that in view of (13.61) . as long as Vn and Y are bounded in time, we have t (f) — 

I) « fV2 ^ 0 as f —)■ 0. By the continuity in time of Vo(f) and ZQ^t), it follows that there exists 


Ti = Ti{M) > 0 


such that 


sup Vo{t) < SCqM 
te[o,Ti] 

(3.9) 

sup Zo{t) < 

tG[0,Ti] ^ 

(3.10) 

This is a time of local existence in for Vu and a. 

At this stage, we assume that T obeys 


T <Ti 

(3.11) 


and define 


Bm — Vm + Z^n — SUp {Ymit) + Z^nit)) 
t£[0,T] 

for all m > 0. By (I3.9I) - (I3.10I) we have 

Bq < 3CqM + —. 

Adding (13.41) and (13.81) we arrive at 

Bm < + CirV2(i + Bo + T^/^Bq + TBl)Bm 


+ CirV2(i + rV2) ^ / .]B,Bm-, 

0<j<m 


m 
j k 


Bj Bf^ B^a—j—k 


0<\{j,k)\<m, 

for all m > 1, for some positive constant (7i > 1. In view of (13.121) we may take 

0 < T = T{Bo) = T{M) < Ti 


sufficiently small, such that 

+ Bo + T^/^Bo + TBl) < 

We thus obtain from (13.131 ) and (13.141 ) that 


Bm < 2Cinm + 2CiT^/^{l + rV2) ^ 

0<j<m 



Bj B^a—j 



Bj Bf^ B^n—j—k 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


for all m > 1. 
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Finally, denote 


||(Vi;,F-/)||5,.,r= 




m>0 


m 


Is 


Multiplying (13.151) by 6"^ml noting that since s > 1 we have (™) ^ < 1 and * < 1, and 

recalling the initial datum assumption (13.11) . we arrive at 

_ p rj jD rm—j 

\\{Vv,Y-I)\\ s,s,t<2CiM + 2CiT^/\1 + T^/^)Y, E ' 


+ 2Cir3/2 Y 


j!® (m — j)!^ 

m>0 0<j<m ' ' 

j!® A:!® (m — j — A:)!® 


m>0 0<\{j,k)\<m 

< 2CiM + 2CiT^/^{l + T^/^)\\{Vv,Y - 
+ 2CiT3/2||(Vn,y-/)||L^. 


(3.16) 


Here we used the discrete Young inequality G In order to conclude the proof, we note 

that the initial values are Vno obeying (13.11) . and Yq = I. Thus, at T = 0 we have 


||(Vi;,y-/)|U,.,o<M, 

and in view of (13.161) . if T is taken sufficiently small so that 

8C,2j.i/2(i 32CfT^/^M‘^ < 


(3.17) 

(3.18) 

(3.19) 


we arrive at 

\\{Vv,Y -I)\\s,s,t<^CiM. 

In summary, we have proven that there exists T = T{M) > 0, given by (13.1 II) . (13.141) . and 
(13.181 ). such that 


^ sup 

^>oIG[0 ,T] 


^T^v{t)\\Hr + 


H/2 


ml 


Is 


<CY\\dT^^o\\Hr 

m>0 


ml 


Is 


= CM 


(3.20) 


for some constant C > 0. This concludes the proof of the a priori estimates needed to establish 
Theorem 1 1.5 1 


Remark 3.1 (Justification of the a priori estimates). Here we show that by using an approx¬ 
imation argument we may rigorously justify the inequality (13.201) . Assume that the initial datum vq 
is real-analytic (e.g., a mollified approximation of fhe original dafum) and if safisfies fhe inequalify 
dHT]), i.e., 

°° rm 

Y\\9T^^o\\Hr—<M (3.21) 

m!® 

m=0 

for some <5 > 0 and s > 1. Then by IIBBZ761lKVllbll we know fhaf fhe solution is real-analytic on 
[0, Ti), where Ti > 0 (cf. (13.111) 1 is fhe time of existence of fhe solution v in which under 

fhe assumptions of fhe fheorem may be faken independenfly of fhe mollification paramefer, and in 
particular if is infinite when d = 2. Thus Bm{t) < oo for all t G [0, Ti) and all m > 0. 
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Let mo > 0 be an arbitrary integer, and define Bm = Bm for m € {0, 1 , , mo} and Bm = 0 
for m G {mo + 1, mo + 2,...}. Similarly, denote by th® same type of truncation corresponding 
to Qm, for all integers m > 0. Then B^, and satisfy the same recursion relation (13.131) . i.e.. 


Bm < CiQm + + Bo + T^/^Bo + TbI)B 


0<j<m 

m 

J k 


+ Y 


m\- 

^ j^m—j 


^ j ^ k ^m—j—k 


for all m > 0. Denote 


Smoit) = Y 

m=0 


Bm(t)d- 


m 


|S 


Bm{t)d” 


m=0 


m 


\S 


Note that Smo is a continuous function of time and 


Smo{0)<M 

Following the derivation in (13.161) . we then obtain 

Smoit) < 2CiM + 2CiTV2(i + + 2CiT^/^sl^ 

for all t > 0. By (13.221) and the continuity of Smo{t)-> we get 


(3.22) 


Smo{T)<^CiM 


(3.23) 


provided that T < Ti is chosen to obey (13.111) . (13.141) . and (13.181) . The bound (13.231) may be 
rewritten as 

Y ^21^6^ < ACiM 
m!® 

m=0 


for all i G [0, T], with T as above. Finally, since mo > 0 is arbitrary, from the monotone conver¬ 
gence theorem we obtain 


Bm{t) ^m 

2-^ m!® 


< 4C7iM 


for all t G [0, T]. Passing to zero in the mollification approximation completes the proof. 


4. Local in time persistence of the Lagrangian Gevrey radius 

In this section we prove Theorem 1 1.1 1 For simplicity of the presentation, we give here the proof 
for d = 3. Fix s > 1 and (5 > 0 so that Vug G Gg^s with norm M, that is, the quantity 

nm:= Y 

|Q|=7n 


obeys 


Y < M 

m!* 


m>0 


(4.1) 
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Fix T > 0, to be chosen later sufficiently small in terms of M, s, and 5. Similarly to the previous 
section for m > 0 define 


VmiT)= sup l|5“Vu(f)||H’-, 

*6[0,r] 1^1^^ 

(4.2) 

Z^{T)= sup Y. ll'9“(E(t)-/)||LU-. 

^6[0,T] \a\=m 

(4.3) 


In order fo esfimafe Vu and ifs derivafives, we use fhe fhree-dimensional curl-div sysfem (11.141) 
and (11.151) fo write 


(curlu)”* = emlkdlV^ = + ^UkiSim - Y^)dlV^ + £rajk{5jl - Y-)div'" 

- EijkiSim - YniSjl - Yj)div^ 

divu = {6ik - Y^)dkv\ 


(4.4) 

(4.5) 


From (I4.5I) - (I4.4I) we conclude fhaf for a G Nq we have 

||5“Vi;||h- < CWd^u^Wnr + C\\d^{£ijk{6im - Yr){6ji - Yj)div>^)\\Hr 

+ C\\d^{e,njk{6ji - Yj)div^)\\Hr + C\\^‘^{eijk{S^m - Yr)djv’^)\\Hr 
+ C\\d^{{6,k-Y^)dkV^)\\Hr. 


Summing fhe above inequality over all multi-indices wifh |q;| = m and faking a supremum over 
t G [0, T] we arrive af 

Vm < CQm + CTZmZoVo + CTZiVm + ZoV^ + CrY^ZmVo 

^ ^ M {t-Y‘^\\df^(Y - I)\\Hr\\d^-^Vv\\Hr) 

0<j<m\a\=m,\l3\=j,l3<a iG[0,T] 

-gt E E G;) 

te[o,r] ^ ^ 

< CQm + CTZ^ZoVo + CTZiVm + CTY^ZoV^ + CTY^Z^Vo 

+ CT^/^ ^ h\z,V^_,+CT (Z)zjZkVm-j-k, (4.6) 


for all m > 1. In (14.61) we have used fhaf if {aa}, {ba}, {ca} are non-negative mulfi-indexed 
sequences, fhen 


E 

|a|=m,|/ 3 |=i,/ 3 <o 



< 



(4.7) 
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and 


E 





< 



(4.8) 


These inequalities follow e.g. from BKVllal Lemma 4.2] and BKVllbl Lemma A.l] and the fact 
that < (j^j)- Indeed, for (14.71) (the proof of (14.81) being analogous), we have by using the 
substitution 7 = a — /3 



(4.9) 


Note that when m = 0, the bound (14.61) reads as 


Vb < CoQo + + Zo)Vo (4.10) 


for some constant Cq > 0. 

As in the two-dimensional case, in order to bound Z^ we appeal to the integral formula for 
Y{t) — I, namely (13.61) . We apply 9“ to identity (13.61) . sum over all multi-indices with \a\ = m, 
divide the resulting inequality by and take a supremum over t G [0, T]. By appealing to (14.71) 
and (14.81) . similarly to (13.81) we obtain 

Zm < CT'/'‘(TZSV„, + TZ„ZoVo + T'/^ZoV„, + + v„) 

+ cr’/= ^ (™)z,z^K„_,_t + crX;‘(’")z,.y„^ (4,ii) 

0<|(i,fc)|<m ^ j=l 

when m > 1 , and 

Zq < CqT^/^{1 + T^/^Zq^Vq (4.12) 


for m = 0 . 

Once the recursive bounds (I4.6I) - (I4.10I) and (14.11I) - (I4.121) have been established, we combine 
them with the initial datum assumption (14.11) . and as in Section [3] obtain that there exists T = 
T{M) > 0 such that 


sup (\\d‘^Vv{t)\\Hr + 

q,>0*S[0,T] \ 


\\d-{Y{t)-I)\\Hr 

tll2 


< 


C^Wd'^VvoWw 


5 l«l 


o>0 


a ! 


Is 


= M, 


(5l“l 

a|!® 


for some constant C > 0. This concludes the proof of the a priori estimates needed to establish 
Theorem 1 1.1 1 
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5. Example of Eulerian ill-posedness in the analytic class Gi 5 

In this section we prove Theorem 11.31 The idea of the proof is similar to the example given 
earlier in Remark 11.21 but addresses the fact that functions whose holomorphic extension have a 
simple pole at ±i(5 do not lie in Gi^s (a fact encoded in the sum over m, as opposed to a supremum 
over m, defining our real-analytic norm, cf. (11.171) and (15.11) 1. To address this issue we integrate such 
a real-valued function four times, so that the holomorphic extension to the strip of radius 6 (where 
(5 = 1) around the real-axis is also a function up to the boundary of this strip (cf. (15.51) ). The 
proof then proceeds by cutting off in a Gaussian way at infinity (cf. (15.91 )). which is compatible with 
real-analyticity, and then periodize the resulting function so that we are dealing with a finite energy 
function (cf. (15.131 )). Verifying that the resulting function ip yields the necessary counterexample to 
prove the theorem follows then from a direct but slightly technical calculation. 

Let /, g be two 27 r-periodic functions. Recall (cf. 11DM871IBT101 ) that the function defined by 

u{xi,X2,X3, t) = (/(X 2 ), 0, g{xi - tf{x2))) 

is an exact solution of the Euler equations posed on T^, where T = [—vr, tt] with the initial datum 

Uo{xi,X2,X3) = {f{x 2 ), 0 ,g{xi)). 

Also, for a 27r-periodic function ip and for (5 > 0 by definition we have that 

( \ a™, 

= E E (5-1) 

7n=0 J 

Note that C C°(T3) in view of the Sobolev embedding. Without loss of generality we fix 

(5 = 1 fhroughout this section. 

We start with a few considerations on the real line M. For a function F € L^(M) we normalize 
the Fourier transform as 

F{0 = ^ / 

v2vr 7 r 

Consider the two decaying real-analytic functions 

hi{x) = \ ^ „ 

^ ^ V vrl + x^ 

and 

These functions have explicit Fourier transforms that are given by 

hi{i) = exp(-l^l) 

and 

^2(0 = exp(-|^| 2 ). 

Define 

h{x) = hi{x) - (1 - + ^(-A) - ^(-A)3/2^ h2{x) 

In view of the above formulae we have that 

m = exp(-ICI) - (1 - ICI + - ^1^1") exp(-leP). 


(5.2) 

(5.3) 

(5.4) 
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Note that 


and 


m = ^^+om^) as lei^o 


24 


HO = exp(-ICI) + O exp - 




as 1^1 —oo. 


Lastly, we define 


so that 


px pxi pX 2 rx 3 

H{x)= / / / / hH0dxAdx^dx2 

Jo Jo Jo Jo 


dx^ 


dx‘^ 


H{x) = h{x). 


By taking the Fourier transform of the above equation we arrive at 


H(n = M = J_ 

{iO^ 


exp(-iei) - ( 1 - lei +-|C|2 --len exp(-ieP) 


Clearly, 


sup |Ff(OI + sup 


ex 


<Co 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


l?l<i l«l>i 

for some constant Cq > 0. The function H however is not in H since it grows as |x| —)• oo, and the 
above computations are formal. To fix this issue, we set 

„2' 

(5.9) 


$(x) = exp (^-y j H{x). 

This function is smooth, and decays as \x\ oo. Moreover, in view of (15.71) and using the explicit 
Fourier transform of the Gaussian, we have 

ii-ri?' 


I>(0 = / exp 
7 r 

/ 


= ; exp 
7 r 

We claim that 


2 

1 


H{r])dr] 
|?7|4 


exp(-|r?|) - ( 1 - |ry| + - '^\r]f ) exp(-|? 7 | 2 ) ) dr]. 


sup |$(0I + sup (lCI'^exp(l^l)ld(Ol) < Cl 

ISI<i ISI>i ^ ^ 

for some universal constant Ci > 0. In order to check whether (15.101) holds, we write 


(5.10) 


|Crexp(|C|)l>(0 = - [ 

Jr 


exp -^- j exp(|,el - hl)l?l 

1 - (1 - |r?| + ||r?|2 - I|r^|3)exp(-|r?p + 


\v\' 


-dr], 


decompose the above integral in the regions 
1 


W<i}, UslolSISp'*}, {|£P'‘‘<H<|{|}. {|i(l>l«|} 


and use both the decay resulting from the Gaussian factor and the decay coming from (E 
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A useful observation that shall be needed below is that we have 


||C|'=exp(-|e|)||^ 2 (K) = —^ 


which by Stirling’s estimate 

(27r)l/2„'^+l/2e-n < ^! < g^n+l/2g-n^ ^ ^ 

yields 

1 


^|||e|'^exp(-|fc|)||^ 2 (R) < 


(5.11) 


(5.12) 


Now, we proceed to construct a periodic function with a finite Gi i norm. First, we build a 
27r-periodic function ip by using the Poisson summation applied to the function <h. More precisely, 
let 


p{x) = $(x — 2rmr). 

m=—oo 

Clearly p is periodic, and its Fourier series coefficients obey 


(5.13) 


(5.14) 


for all A; G Z. 

Therefore, using estimates (15.101) and (15.121) . with the Poisson summation formula, we have that 

— in - 

n! 


Gi.i = 


n>0 


dx' 




7/2 (T) 


< 


< 


C\\p\\h^(t:) + C ^ 


n>5 




dx' 




+ 


L2(T) 


^n+2 


L2(T), 


1 

n! 




n>b 


ni 


< cibii„.,T)+cy; (iii«r<i>(oiiL=(isi>i) + iiier+"i'«)iii=(Ki>i)) ^ 


n>5 


and thus 


Gi,i < C\\p\\h6(j) 


+ Cyy(|||{r-‘exp{-|J|)||y,|4l>i, + |||jr-2exp(-|{|)||t=,|4l>i, 

n>5 ^ 

^ C'||y?||^6(Tr) + C ^ 


1 

n! 


n>5 


(n — 4)! ^ (re — 2)! ^ 1 


< 


c'IIv^II776(t) + C'^ — 


n>5 


(re — 4)1/^ (re —2)^/'^/ re! 

1 


< C^p < oo. 


(5.15) 
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Note that ||(/?||gi ^ = oo for any analyticity radius 5 > 1, since 

(n-2)! <5^ 

_^ ry^ 

(re — 2)^/^ re! 

whenever <5 > 1, and the estimate in (15.151) may also be turned into lower bounds. 

Proof of Theorem 1 1.31 Consider 

g{x) = ip{x) (5.16) 

where (p is as given in (15.131) . and define 

/(x) = sinx (5.17) 

Since / is entire, we have that ||/||gi ^ < oo for any 5 > 0. With the definitions of / and g above, 
it follows from (15.151) that 

IkollGi.i < OO. 

Note that in view of the periodicity in xi and X 2 , the functions f{x 2 ) and g{xi) have finite energy 
(i.e., 77^(T^) becomes 77^(T), up to a multiplicative constant), and the multi-index summation in 
dSD becomes a simple sum over re > 0. Thus (11.191) is established. 

In order to establish (11.201) . we assume, for the sake of obtaining a contradiction, that for some 
t G (0,1/10] we have ||re(f)||Gi i < oo. We fix fhis value of f G (0,1/10] fhroughouf fhis proof. 
Consider fhe function 

ip{xi,X 2 ) := d^^U 3 {xi,X 2 ,X 3 ,t) = g”'{xi - tf{x 2 )). (5.18) 

The inequalify jjre(f)jjGi i < oo implies 

a>0 Ml'/ 

If follows thaf for any R G (0,1), fhe joint in (xi, X 2 ) power series of 1 / at the origin 

l/(xi,X2)= ^ am,nX'^X2 (5.19) 

m,n>0 

converges absolutely in the closed square of side length R at the origin 

Cr = {{xi,X2 ): Ixi] < R, \X2\ < 7?} 

and defines a real-analytic function of fwo variables in fhis square. Thus, we may consider fhe 
complex exfension 

i^{zi,Z2) =ll^{xi+iyi,X2+iy2) = ^ am,nZTz2 

m,n>0 

which converges absolutely when \zi\ < R and \z 2 \ < R- Fix 

Rt = l-j (5.20) 

which clearly belongs fo (0,1), and is thus an allowable choice for R. Also, fix 

xi = 0 and ^;2 = 0 + flog2. 

Since t G (0,1/10), we have \z 2 \ = log 2 < Rt, so fhaf by fhe above consideration, 

lim j'i/(iy 2 ,ilog2)j < oo. 

?/2—/ — 


E 

n>5 


(5.21) 
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In order to complete the proof by contradiction, we shall next show that in fact (15.211) is false, and 
in fact we have 


lim _ |' 0 (iy 2 ,ilog 2 )| = oo. (5.22) 

y2^-Rt 

The remainder of this proof is devoted to establishing (15.221) . 

First observe that sin(i log 2) = 3i/4, and thus 

V^(z2/2,ilog2) = ip'"{i{y2 - 3f/4)). 

Next, note that by the definition of Rt, (15.201) . we have 

3f , 

2/2 J 2/2 —Rt ■ 

Thus, proving (15.221) amounts to showing that 

lim \ip"'{iy)\ = oo (5.23) 

y^-l+ 

which is what we establish below. In view of (15.91 ). (15.131) . and the Leibniz rule, we have that 

^"'{z - 2rmr) (5.24) 

mez\{0} 

and 

4>"'(z) = exp - 3zH"{z) + 3{z^ - l)H'{z) + z{3 - z^)H{z)) (5.25) 

for any complex number z with \z\ < 1. Next, note that by (15.41) and (15.61) . we have 

H^^^\z) = hiz)+£oiz) 


where 

£oiz) = (l - (-A)V2 + ^(_A) - ^(-A)3/2^ h 2 {z) (5.26) 

is an entire function (since its Fourier coefficients are given by a polynomial times a decaying 
Gaussian). Moreover, letting £^ 1 ( 2 ;) = £o(wi)dwi, £ 2 (z) = £o(w 2 )diV 2 dwi, £ 3 (z) = 

fofo^^fo^^ ^o(w 3 )di 03 dw 2 dwi, and £^ 4 ( 2 :) = £o{w 4 )dw 4 dws dw 2 dwi, we immedi¬ 

ately obtain that 

£{z) = £i{z) - 3 z£2{z) + 3 ( 2 ^ - 1 )^ 3 (- 2 ) + ^(3 - z'^)£ 4 {z) (5.27) 


is also an entire function. On the other hand, we may explicitly compute the integrals of h as 


Riiz) = j hi{wi)dwi = 


arctan 2: 


n W\ 

hi{w2)dw2 dwi = 


2 arctan z — - log(l + 2^) 


n wi rW2 

/ hi{w3)dw3 dw2 dwi = -1/ — (2 + (2^ — 1) arctan 2 — 2log(l + 2^)) 
7o 2 V TT 

n W\ rW2 /*U73 

/ / /ll {W4)dw4 dW 3 dW 2 dwi 

Jo Jo 

= — \ — (52^ + 22(2^ — 3 ) arctan 2 — ( 32 ^ — 1 ) log(l + 2^)) 

12 V TT ^ ' 
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which implies that 

niz) := ni{z) - 3 z'H2{z) + 3(^2 - l)'H 3 iz) + z{3 - z^'j'H^iz) 

+ 2(15 — 45z^ + — z®) arctan z + 2;(39 — 28z‘^ + 3z^) log(l + . (5.28) 

In summary, with the definition of S in (15.271) and of Ti in (15.281) . we have that 

4>"'(z) = exp (^(^) + ^(^))- (5-29) 

Letting 2 = iy, and using that arctan(fy) = iarctany, we arrive at 

$"'(iy) = exp {T-Liiy) + £{iy)). (5.30) 

Since £ is an entire function, we have that supj^g[_i gj \£{iy)\ < C < 00 . Writing 

niiy) = fW-y^arctanhy + -^W-y(-18 - 33y'^ - 5y^) 

+ Y (39 + 28y^ + 3 /) (2 arctanh y + y log(l - y‘^)J 

+ ^a/-(- 24 + lly^ + 12y^ + y®) arctanh y (5.31) 

6 V vr 

and observing limy_^_i+ (2 arctanh y + y log(l — y^)) = — log 4 and limj^_^_i+ (y+1) arctanh y = 
0, we arrive at 


lim \'H{iy)\ = 00 (5.32) 

y^-l+ 

since arctanh has a logarithmic singularity at y = —1. Combined with the above, it follows from 
(15.321 ) that 


lim |4>"'(fy)| = 00 (5.33) 

y^-l+ 


which in turn shall imply that (15.231) holds. 

Indeed, the only remaining part of the proof is to show that 


lim 

y^-l+ 


\^'"{iy — 2m7r)| < 00 . 


mSZ\{0} 

The above holds since for each m / 0 we have that 

|7f(fy — 2m7r)| + \£{iy — 2m7r)| < P{m) 
uniformly for |y| G [1/2,1], where P is a polynomial, and since 


exp 


{iy — 2rmry 


< exp 


-2m^7r^ 


which makes the sum over m / 0 finite. In order to obtain the first bound, we use (15.311) and the 
formula 

arctan z = -f(log(l — iz) — log(l + iz)') 

where the complex domains of the above logarithms are cut on [0, 00 ) and (— 00 ,0] respectively. □ 
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